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Abstract: 

We present an improved calculation of -B — light vector form factors from light-cone sum 
rules, including one-loop radiative corrections to twist-2 and twist-3 contributions, and lead- 
ing order twist-4 corrections. The total theoretical uncertainty of our results at zero mo- 
mentum transfer is typically 10% and can be improved, at least in part, by reducing the 
uncertainty of hadronic input parameters. We present our results in a way which details 
the dependence of the form factors on these parameters and facilitates the incorporation of 
future updates of their values from e.g. lattice calculations. We also give simple and easy- 
to-inplement parametrizations of the g^-dependence of the form factors which are valid in 
the full kinematical regime of q^. 
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1 Introduction 



This paper aims to give a new and more precise determination of the decay form factors 
of Bd^u,s mesons into hght vector mesons, i.e. p, u, K* and 0; it is a continuation of of our 
previous study of B decays into pseudoscalar mesons [1]. The calculation uses the method 
of QCD sum rules on the light-cone, which in the past has been rather successfully applied 
to various problems in heavy-meson physics, cf. Refs. [2];^ an outline of the method will be 
given below. Our calculation improves on our previous paper [4] by 

• including B ^ ui form factors; 

• including radiative corrections to 2-particle twist-3 contributions to one-loop accuracy; 

• a new parametrization of the dominant hadronic contributions (twist-2 distribution 
amplitudes) ; 

• detailing the dependence of form factors on distribution amplitudes; 

• a new parametrization of the dependence of the form factors on momentum transfer; 

• a careful analysis of the theoretical uncertainties. 

Like in Ref . [1] , the motivation for this study is twofold and relates to the overall aim of B 
physics to provide precision determinations of quark flavor mixing parameters in the Stan- 
dard Model. Quark flavor mixing is governed by the unitary CKM matrix which depends 
on four parameters: three angles and one phase. The constraints from unitarity can be 
visualized by the so-called unitarity triangles (UT); the one that is relevant for B physics 
is under intense experimental study. The over-determination of the sides and angles of this 
triangle from a multitude of processes will answer the question whether there is new physics 
in flavor-changing processes and where it manifests itself. One of the sides of the UT is 
given by the ratio of CKM matrix elements iKib/Kfel- Wcb\ is known to about 2% accuracy 
from both inclusive and exclusive h — > civ transitions [5], whereas the present error on \Vuh\ 
is much larger and around 15%. Its reduction requires an improvement of experimental 
statistics, which is under way at the B factories BaBar and Belle, but also and in particular 
an improvement of the theoretical prediction for associated semileptonic spectra and decay 
rates. This is one motivation for our study of the B ^ p semileptonic decay form factors Ai, 
A2, V, which, in conjunction with alternative calculations, hopefully from lattice, will help 
to reduce the uncertainty from exclusive semileptonic determinations of \ Vub\- Secondly, form 
factors of general B hght meson transitions are also needed as ingredients in the analysis 
of nonleptonic two-body B decays, e.g. B pn, in the framework of QCD factorization 
[6], again with the objective to extract CKM parameters. One issue calling for particular 
attention in this context is the effect of SU(3) breaking, which enters both the form factors 
and the K* and meson distribution amplitudes figuring in the factorization analysis. We 
would like to point out that the implementation of SU(3) breaking in the light-cone sum 
rules approach to form factors is precisely the same as in QCD factorization and is encoded 
in the difference between p, cu, K* and 4> distribution amplitudes, so that the use of form 

^See also Ref. [3] for reviews. 
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factors calculated from light-cone sum rules together with the corresponding meson distri- 
bution amplitudes in factorization formulas allows a unified and controlled approach to the 
assessment of SU(3) breaking effects in nonleptonic B decays. 

As we shall detail below, QCD sum rules on the light-cone allow the calculation of 
form factors in a kinematic regime where the final state meson has large energy in the 
rest-system of the decaying B, E ^ Aqcd- The physics underlying B decays into light 
mesons at large momentum transfer can be understood qualitatively in the framework of 
hard exclusive QCD processes, pioneered by Brodsky and Lepage et al. [7]. The hard scale 
in B decays is and one can show that to leading order in l/rub the decay is described 
by two different parton configurations: one where all quarks have large momenta and the 
momentum transfer happens via the exchange of a hard gluon, the so-called hard-gluon 
exchange, and a second one where one quark is soft and does interact with the other partons 
only via soft-gluon exchange, the so-called soft or Feynman-mechanism. The consistent 
treatment of both effects in a framework based on factorization, i.e. the clean separation 
of perturbatively calculable hard contributions from nonperturbativc "wave functions", is 
highly nontrivial and has spurred the development of SCET, an effective field theory which 
aims to separate the two relevant large mass scales nib and ^/rn^JiQCD in a systematic way 
[8] . In this approach form factors can indeed be split into a calculable factorizable part which 
roughly corresponds to the hard-gluon exchange contributions, and a nonfactorizable one, 
which includes the soft contributions and cannot be calculated within the SCET framework 
[9, 10]. Predictions obtained in this approach then typically aim to eliminate the soft part 
and take the form of relations between two or more form factors whose difference is expressed 
in terms of factorizable contributions. 

The above discTission highlights the need for a calculational method that allows numerical 
predictions while treating both hard and soft contributions on the same footing. It is precisely 
QCD sum rules on the light-cone (LCSRs) that accomplish this task. LCSRs can be viewed 
as an extension of the original method of QCD sum rules devised by Shifman, Vainshtein and 
Zakharov (SVZ) [11], which was designed to determine properties of ground-state hadrons at 
zero or low momentum transfer, to the regime of large momentum transfer. QCD sum rules 
combine the concepts of operator product expansion, dispersive representations of correlation 
functions and quark-hadron duality in an ingenious way that allows the calculation of the 
properties of non-excited hadron-states with a very reasonable theoretical uncertainty. In the 
context of weak-decay form factors, the basic quantity is the correlation function of the weak 
current and a current with the quantum numbers of the B meson, evaluated between the 
vacuum and a light meson. For large (negative) virtualities of these currents, the correlation 
function is, in coordinate-space, dominated by distances close to the light-cone and can 
be discussed in the framework of light-cone expansion. In contrast to the short-distance 
expansion employed by conventional QCD sum rules a la SVZ where nonperturbativc effects 
are encoded in vacuum expectation values of local operators with vacuum quantum numbers, 
the condensates, LCSRs rely on the factorization of the underlying correlation function into 
genuinely nonperturbativc and universal hadron distribution amplitudes (DAs) which are 
convoluted with pro cess- dependent amplitudes T. The latter are the analogues of Wilson- 
coefficients in the short-distance expansion and can be calculated in perturbation theory. 
The hght-cone expansion then reads, schematically: 




(1) 



n 
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The sum runs over contributions with increasing twist, labelled by n, which are suppressed 
by increasing powers of, roughly speaking, the virtualities of the involved currents. The 
same correlation function can, on the other hand, be written as a dispersion-relation, in the 
virtuality of the current coupling to the B meson. Equating dispersion-representation and 
the light-cone expansion, and separating the B meson contribution from that of higher one- 
and multi-particle states using quark-hadron duality, one obtains a relation for the form 
factor describing the decay B — > hght meson. 

A crucial question is the accuracy of light-cone sum rules. Like with most other methods, 
there are uncertainties induced by external parameters like quark masses and hadronic pa- 
rameters and intrinsic uncertainties induced by the approximations inherent in the method. 
As we shall discuss in Sec. 4, the total theoretical uncertainty for the form factors at = 
is presently around 10%, including a 7% irreducible systematic uncertainty. 

Our paper is organized as follows: in Sec. 2 we define all relevant quantities, in particular 
the meson distribution amplitudes. In Sec. 3 we outline the calculation. In Sec. 4 we 
derive the sum rules and present numerical results which arc summarized in Tabs. 7 and 
8. Section 5 contains a summary and conclusions. Detailed expressions for distribution 
amplitudes, a break-down of the light-cone sum rule results into different contributions and 
explicit formulas for the contributions of 3-particle states are given in the appendices. 



2 Definitions 

B ^ V transitions, where V stands for the vector mesons p, oo, K* and 0, can manifest 
themselves as semilcptonic decays B Viug or rare flavor- changing neutral current (FCNC) 
penguin-induced decays B V'j and B —>■ Vi'^i~. All these decays are described by a total 
of seven independent form factors which usually are defined as {q — Pb — p) 



Cv{V(p)\q-fi,(l - -f5)b\B(pB)) = -iel(mB + mv)Ai(q ) + i(pB + p)i,(e*q) 



TTlB + my 



+ zq,{e*q) ^ {A,{q') - Ao{q')) + e,.,,e*^p%p'^ (3) 



with A,{q') = + A,{q') - VlE^ A,{q^) and A,{Q) = ^3(0); (3) 

Cv{V{p)\qa^,q''{l + ^^)h\B{pB)) = ie^,pae*''p'BP'"^T^{q^) 

+ r2(g') [eliml - ml) - {e*q) {pb+p)^.} + n{q^){e*q) - , , {Pb+p)A (4) 

withri(O) =^2(0). (5) 

Ao is also the form factor of the pscudoscalar current: 

cviVld^A^'lB) = cvirrib + mg){V\qi^5b\B) = 2mv{e*q)Ao{q^). (6) 

q in the above formulas stands for u, d and s in (2) and (6) and d, s in (4); the actual 
assignment of different decay channels to underlying b ^ q transitions is made explicit in 
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Tab. 1. In our calculation, we assume isospin symmetry throughout, which implies that 
there are five different sets of form factors: Bq ^ p, Bq ^ lu, Bq K*, Bg ^ K* and 

— > (with q — u,d). The factor cy accounts for the fiavor content of particles: cy — 
for uj and cy = 1 otherwise.^ 

The currents in (2) and (4) contain both vector and axialvector components. V and Ti 
correspond to the vector components of the currents, and, as the B meson is a pseudoscalar, 
to the axialvector components of the matrix elements. Ai^2 clearly correspond to the ax- 
ialvector component of the V — A current; the term in A2, — Aq arises as the contraction 
of (2) with —iq^^ must agree with (6). As for the penguin current, T2 3 correspond to the 
axialvector components of the current; there is no analogon to Aq, as the current vanishes 
upon contraction with q^. As we shall see in Sec. 4, for analysing the dependence of each 
form factor on (f, it is best to choose ^0,1,2 as independent form factors for the A current, 
and define A3 by (3), but for the penguin current it will turn out more appropriate to choose 
a different set of independent form factors: Ti, T2 and Ta with 

Cv{V{p)\qa^,q\l + l^)h\B{pB)) = ^e^.pae»<^ 2T^{q^) 

+ e;(m| - ml)T2{q^) - {ps + p),{e*q)Ts{q') + q,{e* q)Ts{q') (7) 

and T3 defined as 

Ts{q') = '^^^^{Uq')-U^^^^ (8) 

As the actual calculation is done using an off-shell momentum ps with p^ ^ it is 
crucial to avoid any ambiguity in the interpretation of scalar products like 2pq = p%—q^—p^ 7^ 
m?^ — q^ — my that occur at intermediate steps of the calculation. This is particularly relevant 
for the penguin form factors which arc defined in terms of a matrix clement over the tensor 
current which is contracted with the physical momentum q^ . The problem can be avoided 
by extracting Tj and T3 from sum rules for a matrix element with no contractions: 

{V{p)\qa^,u-i^h\B{pB)) = A{q^) {e;(pB - (pb - B{q^) {e*^q, - q^el] 

- 2C(g2) — 2" {Pi^^l^ - %P^) ■ (9) 

A, B and C are related to Tj and T3 defined in (4) and (7) as 

Ti = A, T2 = A- 2 ^ B, Ts = B + C, T3 = A + ^ ^ C, (10) 

— my m^ — ruy 

which implies 

Ti(0)= 7^2(0) =73(0). (11) 

Relevant for semileptonic decays are, in the limit of vanishing lepton mass, the form 
factors A12 and V with g^, the invariant mass of the lepton-pair, in the range < < 
{mB — mvY- B — > depends on Ti(0), whereas B — > Vt^l~ depends on all seven form 

To be precise, cy is V2 for in 6 ^ M and for lo, and — \/2 for in 6 ^ d, with the flavor wave 
functions p° ~ {uu — dd) / and uj ~ {uu + dd) j \/2. We assume that is a pure ss state. 
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p+ p^.oj p- K*+ K*%ds) K*- K*%sd) (j) 

— b ^ u b ^ d — — b ^ s — — 

Bd b ^ u b ^ d — — — — b s — 

Bs — — — b ^ u b ^ d — — 6— >s 



Table 1: Allowed decay channels in terms of underlying quark transitions. We assume 
isospin-symmetry and hence have five different sets of form factors: Bg ^ p, Bg ^ cu, 
Bg K*, B, K* and B.^cf) (with q^u,d). 



factors (see Ref. [12] for an explicit formula). The motivation for studying B — > pii/i and 
B uolvf is to extract information on the CKM matrix element \Vuh\-, whereas the FCNC 
transitions B {K* , p,uj)'~f and B {K* , p,uj)i'^£^ serve to constrain new physics or, in 
the absence thereof, the ratio |Vts/V(d| [13], which would complement the determination of 
\Vt.s/Vtd\ from S-mixing. 

In the LCSR approach the form factors are extracted from the correlation function of the 
relevant weak current Jw, i-c. cither the semileptonic V — A current or the penguin current 
of (8), and an interpolating field for the B meson, in the presence of the vector meson: 

^(g^p|) = i J d'xe"^^{V{p)\TMx)jlmO), (12) 

with jh — 'mhq'ij5b, q' e {u, d, s}. For virtualities 

ml-pl> ©(AqcDrnft), ml-q^> 0{KqcT>mj,), (13) 

the correlation function (12) is dominated by light-hke distances and therefore accessible to 
an expansion around the light-cone. The above conditions can be understood by demanding 
that the exponential factor in (12) vary only slowly. The light-cone expansion is performed 
by integrating out the transverse and "minus" degrees of freedom and leaving only the 
longitudinal momenta of the partons as relevant degrees of freedom. The integration over 
transverse momenta is done up to a cutoff, /xir, all momenta below which are included 
in a so-called hadron distribution amplitude (DA) 0, whereas larger transverse momenta 
arc calculated in perturbation theory. The correlation function is hence decomposed, or 
factorizcd, into pcrturbative contributions T and nonperturbative contributions 0, which 
both depend on the longitudinal parton momenta and the factorization scale /^ir. If the 
vector meson is an effective quark-antiquark bound state, as is the case to leading order in 
the light-cone expansion, one can write the corresponding longitudinal momenta as up and 
(1 — ii)p. where p is the momentum of the meson and u a number between and 1. The 
schematic relation (1) can then be written in more explicit form as 

n<l\pl)=Y. f duT^-\u,q\pl,p,^)ct>^-\u,pi^). (14) 

As r itself is independent of the arbitrary scale /xjr, the scale-dependence of T^") and 0^"^ 
must cancel each other. ^ If 0^"^ describes the meson in a 2-parton state, it is called a 2-particle 

^If there is more than one contribution of a given twist, they will mix under a change of the factorization 
scale /UiR and it is only in the sum of all such contributions that the residual /xir dependence cancels. 
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DA, if it describes a 3-parton, i.e. quark-antiquark-gluon state, it is called 3-particle DA. In 
the latter case the integration over u gets replaced by an integration over two independent 
momentum fractions, say ai and 0:2 ■ Eq. (14) is called a "coUinear" factorization formula, 
as the momenta of the partons in the meson are coUinear with its momentum. Any such 
factorisation formula requires verification by explicit calculation; we will come back to that 
issue in the next section. 

Let us now define the distribution amplitudes to be used in this paper. All definitions 
and formulas are well-known and can be found in Ref. [14]. In general, the distribution 
amplitudes we are interested in are related to nonlocal matrix elements of type^ 

{o\q2{o)r[o,xMx)\v{p)) or (o|^2(o)[o, HrG^.MA72K 

X is light-like or close to light- like and the light-cone expansion is an expansion in x^; ^; is a 

number between and 1 and F a combination of Dirac matrices. The expressions [0,x] etc. 
denote Wilson lines that render the matrix elements, and hence the DAs, gauge-invariant. 
One usually works in the convenient Fock-Schwinger gauge x^A1j^{x)X"' /2 = 0, where all 
Wilson lines are just 1; we will suppress them from now on. 

The DAs are formally ordered by twist, i.e. the difference between spin and dimension 
of the corresponding operators. In this paper we take into account 2- and 3-particle DAs 
of twist-2, 3 and 4. The classification scheme of vector meson DAs is more involved than 
that for pseudoscalars; it has been studied in detail in Ref. [14]. One important point 
is the distinction between chiral-even and chiral-odd operators, i.e. those with an odd or 
even number of 7^-matrices. In the limit of massless quarks the DAs associated with these 
operators form two completely separate classes that do not mix under a change of /xir. 
One more important parameter is the polarisation-state of the meson, longitudinal (||) or 
transverse (±), which helps to classify twist-2 and 3 DAs. Up to twist-4 accuracy, we have 
the following decomposition of chiral-even 2-particle DAs [14]: 

P, I du e--^^ [011 (u) + ^ A|| (u) + 0{x')\ 
-^^^.^m^^'c^ne--^^ [g,{u) + <P\\{u)-2g^l\u) + 0{x^)) 1, (15) 



(0|g2(0)7M75gi(^)|^(P, A)) = -- fvmve;^Pe^^^p^zp due-'^^^ (^f V) + 0{x')) , (16) 
and for the chiral-odd ones: 

{el^^P.-ei'^P,) fdu. 

Jo 



*The currents to use for and ui are {uTu =F dTd) / \/2, respectively. 
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(17) 



{OMO)q,{x)\V{P, A)) = I (eW^) m'y £ di^e"^"^^ (/i[f)(i.) + 0{x')) . (18) 

The relevant 3-particle DAs are defined in App. B. 

Note that we distinguish between hght-hke vectors p, z with — — z"^ and the vectors 
P, X with = ruy and x"^ ^ 0; exphcit relations between these vectors are given in App. B. 
The DAs are dimensionless functions of u and describe the probability amplitudes to find 
the vector meson V in a state with minimal number of constituents — quark and antiquark 
— which carry momentum fractions u (quark) and 1 — u (antiquark), respectively. The eight 
DAs (f) — {4>\\,±, g±'°'\ h\f'^\ hs, gs} are normalized as 



^0 



du(/){u) = l. (19) 



The nonlocal operators on the left-hand side are renormalized at scale /x, so that the distri- 
bution amplitudes depend on fi as well. This dependence can be calculated in perturbative 
QCD; we will come back to that point below. 

The vector and tensor decay constants fv and fy featuring in Eqs. (15) and (18) are 
defined as 

(0|g2(0)7M?i(0)|^(P, A)) = fvmyeP, (20) 
(0|g2(0)a^.gi(0)|\/(P, A)) = ^/^'(/i)(e(^)p. - e^P^); (21) 
numerical values are given in Tab. 3. fy depends on the renormalization scale as 

with L = as{Q'^)/as{fJ'^) and /?o = 11 - 2/3n/, Uf being the number of flavors involved. 

The DAs as defined above do actually not all correspond to matrix elements of operators 
with definite twist: ^^^y are of twist-2, /iy*'*-' and fif^'"^ contain a mixture of twist-2 and 
3 contributions and Aj^ y, and a mixture of twist-2, 3 and 4 contributions. Rather 
than as matrix elements of operators with definite twist, the DAs are defined as matrix 
elements of operators built from fields with a fixed spin-projection onto the light-cone. For 
quark fields, the possible spin projections are s = ±1/2 and the corresponding projection 
operators P+ = l/{2pz)'^^ and P_ = l/{2pz)^'^. Fields with fixed spin-projection have a 
definite conformal spin, given by j = 1/2 (s -|- canonical mass dimension), and composite 
operators built from such fields can be expanded in terms of increasing conformal spin.^ The 
expansion of the corresponding DAs, suitably dubbed conformal expansion, is one of the 
primary tools in the analysis of meson DAs, and together with the use of the QCD equations 
of motion it allows one to parametrize the plethora of 2- and 3-particle DAs in terms of 
a manageable number of independent hadronic matrix elements. DAs defined as matrix 
elements of operators with definite twist, on the other hand, do not have a well-defined 
conformal expansion [16], and this is the reason why we prefer the above definitions. In an 

^For a more detailed discussion we refer to the first reference in [14] and to Ref. [15]. 
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admittedly rather sloppy way we will from now on refer to g^^'""^ , /ly'''*'' as twist-3 DAs and to 
/is, gs, A_L^|| as twist-4 DAs. A more detailed discussion of the relations between the different 
DAs is given in App. B; the upshot is that the 18 twist- 2, 3 and 4 DAs we shall take into 
account can be paramatrized, to NLO in the conformal expansion, in terms of 10 hadronic 
matrix elements, most of which give only tiny contributions to the LCSRs for form factors. 

For the leading twist-2 DAs (j)\\^± in particular, the conformal expansion goes in terms of 
Gegenbauer polynomials: 

<P{u, /x^) = 6u{l - ^1 + an{fi^)Cl/\2u - 1) j . (22) 

The first term on the right-hand side, 6u{l — u), is referred to as asymptotic DA; as the 
anomalous dimensions of a„ are positive, (j) approaches the asymptotic DA in the limit //^ — > 
oo. The usefulness of this expansion manifests itself in the fact that, to leading logarithmic 
accuracy, the (nonperturbative) Gegenbauer moments renormalize multiplicatively with 

a„(g^)=L^"/(^^°)a„(//2) (23) 
with L = as{Q'^) / Cis{jJi^)- The anomalous dimensions 7!'^ are given by 

7^ = 8Cf(i*(» + 2)+7e-1) (25) 

with ■?/'(n + 1) = X]fe=il/^ ~ lE- As the contributions from different comformal spin do 
not mix under renormalization, at least to leading logarithmic accuracy, one can construct 
models for DAs by truncating the expansion at a fixed order. Despite the absence of any 
"small parameter" in that expansion, the truncation is justified inasmuch as one is interested 
in physical amplitudes rather than the DA itself. If we write 



amplitude = / du(f){u)T(u), 
Jo 

then, assuming that T is a regular function of u, i.e. with no (endpoint) singularities, the 
highly oscillating behavior of the Gegenbauer polynomials suppresses contributions from 
higher orders in the conformal expansion. Even for a function T with a mild endpoint 
singularity, for instance T — \nu, we find, using the generating function of the Gegenbauer 
polynomials, 

du (j){u)T{u) ^ --ao + y — — — 3a„. 

D n(n + 3) 

n=l ^ ' 

This result indicates that, assuming the a„ fall off in n, which, as we shall see seen in Sec. 4.3, 
is indeed the case, even a truncation after the first few terms should give a reasonable 
approximation to the full amplitude. A more thorough discussion of the convergence of 
the conformal expansion for physical amplitudes can be found in Ref. [17]. The major 
shortcoming of models based on the truncation of the conformal expansion is the fact that 
the information available on the actual values of the a„ (and in particular their analogues 



^0 



8 



in 3-particle DAs) is, to put it mildly, scarce. We therefore use truncated models only for 
DAs whose contribution to the LCSRs is small as is the case for all 3-particle DAs and the 
twist-4 DAs; explicit formulas are given in App. B. All contributions due to or induced by 
twist-2 DAs, on the other hand, are treated as described in Sec. 4.3. 

The major difference between the analysis of LCSRs for B vector meson form factors 
and that of B — > pseudoscalar form factors presented in [1] is probably the identification 
of a suitable parameter by which to order the relative weight of different contributions to 
the sum rules. For B — > pseudoscalar form factors, the standard classification in terms of 
increasing twist proved to be suitable, as the chiral parities of the twist-2 DA and 2-particle 
twist-3 DAs arc different, so that contribution of the latter to the LCSRs is suppressed by a 
factor m\/{mu + md)/mh. In addition, the admixture of twist-2 matrix elements to twist-3 
DAs and of twist-2 and 3 matrix elements to twist-4 DAs is small and moreover vanishes in 
the chiral limit uIt^ — > 0. For vector mesons, the situation is more complex: for instance, 
both the twist-2 DA (f)j_ and the twist-3 DA g^f^ contribute at the same order to the form 
factors A2 and Aq, in the combination 0|| — g^f^ . Naive twist-counting is evidently not very 
appropriate for classifying the relative size of contributions of different DAs to the form 
factors. Instead, we decide to classify the relevance of contributions to the LCSRs not by 
twist, but by a parameter 5 oc mv- The precise definition of 5 depends on the kinematics of 
the process; to leading order in an expansion in l/wb, however, one finds 5hql = mv/rrib. 
The numerical analysis of the LCSRs docs indeed display a clear suppression of terms in 
0{5) and higher, which suggests the following classification of 2-particle DAs: 





■ 0±; 


• 0{6') 






:/i[f'*\/i3,Ax; 


• 0{6^) 


: c/3, A||. 



Wc treat 6 as expansion parameter of the light-cone expansion and shall combine it with the 
perturbative QCD expansion in ag to obtain a second order expression for the correlation 
functions (12); terms in 5^ are dropped. 

3 Calculation of the Correlation Functions 

As we have seen in Sec. 2, LCSRs for form factors are extracted from the correlation function 

of the corresponding weak current with the pseudoscalar current jb = mhq'i'-f^b, evaluated 
between the vacuum and the vector meson. In this section we describe the calculation of 
these correlation functions to second order in ct^ and S. 
The relevant correlation functions are defined as 

i J d^xe'^^ViprnV - A),(x)jlmO) ^ 

= -iPoe; + ir+{e*q) {q + 2p)^ + ir_ {e*q) q^ + Fye/^^e^g^p^ , (26) 
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t j d'xe^''%V{p)\T[qa,.^^h]{x)jl{Qm ^ 

= A{el{2p + q), - el{2p + q)^} - B{e;q, - eM - 2C{e*q){p^q, - q,p,}. (27) 

The definitions of r='= and C differ from those used in Ref. [4] by a factor pq; we shall come 
back to this point below. In this section we describe the calculation of the contributions of 
2-particle DAs to the above correlation functions; those of 3-particle DAs are calculated in 
App. C. 

In light-cone expansion and including only contributions from 2-particle Fock-states of 
the mesons, each of the seven invariants ro,±,y, A,B,C can be written as a convolution 
integral of type 

r^ = / ^<I>1(K)tUK,pI,q') (28) 

with a, b being spinor indices, p^ = my is set to and is the longitudinal momentum of 
the quark in the vector meson V, which is related to the momentum fraction u introduced 
in Sec. 2 by A;+ = up+.^ The above factorisation formula imphes a complete decouphng 
of long-distance QCD effects, encoded in the DA 0^, and short-distance effects calculable 
in perturbation theory, described by T. Factorisation also makes it possible to calculate T 
in a convenient way: if it holds, T must be independent of the specific properties of the 
external hadron state, and one can calculate with a particularly simple state that allows 
a straightforward extraction of the short-distance amphtudes TJ A convenient choice of the 
external state is a free quark-antiquark pair with longitudinal momenta up and up and spins 
s and r, respectively, and DA 



e/^(M = / c^^-e-^'=+^-(gi(«p,5)g2(«P,r)|(gi)„(^)[^,0]fe)6(0)|0) 



z+=0,z±=0 



2Tiuf{up,s)vf{up,r)5{k^ - up+), 



where u and v are the standard fermion spinors. The T amplitudes, to one loop accuracy, 
are then given directly by the diagrams shown in Fig. 1 with external on-shell quarks with 
momenta up and (1 — u)p, respectively. The projection onto a specific Dirac structure is 
done using the general decomposition 

1_1_ l_ 1 

{qi)a{q2)b = I (l)6a(?l?2) " | (hb) ba{qih5q2) + | n) ba{qi^'^ q2) " ^ (7M75)6a(?l 7^75^2) 

+ ^K.)6afe(T^"g2). (29) 

In order to obtain the convolution integrals for vector mesons, one has to replace the struc- 
tures ^iFq'2 in (29) by the appropriate DAs and include factors of e*z, pz and as given in 

^Thc plus-component of a 4- vector k'^ is defined as fc+ = + k^)/y/2, the minus-component as fc_ = 

^This is completely analogous to the calculation of Wilson-coefficients in a local operator product ex- 
pansion, which must be independent of the external states and hence are calculated using any convenient 
state. 



10 







Figure 1: Diagrams for 2-particle correlation functions. V is the weak interaction vertex. The 
hght-quark self-energy diagrams of type c give purely divergent contributions in 1/ eiR— 1 / euy. 



Eqs. (15) to (18). The translation of explicit terms in Zf^ into momentum space is given by 

. d 



—I 



(30) 



as the outgoing comes with a factor ex^{iupz). Terms in l/{pz) can be treated by partial 
integration: 

1 r 

(31) 

— g^f^ , one has 



— Mu) —i dv(j)(v) = —i^(u). 
pz Jo 



There are no surface terms, as for all the relevant structures 4>, e.g. 
$(0) = = ^(1)- A second, approximate way to deal with factors {e*z)/{pz) is based on the 
observation that {e*z) projects onto the longitudinal polarisation state of the vector meson, 
cf. Eq. (A. 3), and that in the ultrarelativistic limit Ey oo the longitudinal polarisation 
vector is approximately coUinear with the meson's momentum: 



.(0) 



mv 



e*z ^1 1 ^ e*q 

pz mv my pq 
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Up to corrections in my, this procedure yields results identical with those from partial 
integration — provided that the corresponding DA is normahzed to 0. That is: 

where the first relation is valid only if (f) is normahzed to 0, i.e. $(1) = 0. This is indeed the 
case for the mixed-twist structure 0|| — g^f^ , but does not apply if only the pure twist-2 DA 
(/)|[ is included, as done in [4]. In this case, unphysical singularities in = appear 
in F-i- and C and have to be factored out. This explains the appearance of additional 
factors l/{pq) = '^/{Pb ~ q'^) in the correlation functions used in [4]. In the calculation 
presented in this paper we use the prescriptions (30) and (31) througout and hence avoid 
unphysical singularities in p^. We have checked that indeed the singularity structure of all 
seven invariants ro,±,y. A, B, C is given by a cut on the real axis for > mf. 

The complete correlation function, including 2-particle DAs to 0(5^) for the vector cur- 
rent V, axialvector current A, tensor current T and scalar current S,^ can now be written 
as 

r(g',p|)- ['du^^,{u;i,')TUu,pl,q';fi') (32) 

C=V,A,T,S 

With = ^ fvmy{r)a, {-P.e} $(«) + e;^^^^)) , (33) 

^1 = -^fvmv{l''l,)abea.x0e*-P'^gt\u)^^^, (34) 

= 4 fvi^^\. {<P, [Mu) - ^ ^^.A. (35) 
+<Pae;/,(.) -\mlelHM ^} , (36) 

= -\ mlf^{<){lUh^\u) (37) 
where $(tt) — J dv (j>\\{v) — g'f\v)J , 

Il{u) = y dv j dw f /i||*^(w) - - (l)±{w) - -h3{w) 



Hsiu) = / dv{h{v)-M^))- 
Jo 

All these three functions F{u) fulfiU F(0) = 1 = F(l). 

^The matrix elements of vector mesons over the pseudoscalar current vanish. 
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Just to give an example, the tree-contribution is given by 

nr = -i [n<i +ur+ mM,a /((? + ^pf - mi) 

with the weak vertex F. 

In order for factorisation to hold, two conditions have to be met: 

(a) the long-distance infrared sensitive parts (IR-singularities) in T have to cancel against 
those in the DAs; 

(b) the convolution integral J du4>'^^"'{u)T'^^"'{u) has to converge. Otherwise factorisation is 
violated by soft end-point singularities. 

In order to check condition (a), we decompose the bare amplitude into finite and divergent 
terms as 

Ultra-violet divergences, which only occur for the penguin-current, are easily subtracted 
using the known renormalisation of the corresponding current: 

T^'^'^u) T^'^'^u) - SZp^^gT^^\u). 

The remaining divergent terms have to cancel against the divergent parts of the bare DA, 

SO that ^ 

(l)'^"{u) is known explicitly for the twist-2 tt DA [7] and coincides with that for 0j|, but to the 
best of our knowledge has not yet been calculated for 0^. Alternatively, one can check the 
cancellation of divergences order by order in the conformal expansion of the DAs, cf. Sec. 2 
and App. C, with^ 

We find that all 1/e terms cancel as required. 

As for condition (b), we also find that all T are regular at the endpoints, so that there 
are no endpoint singularities in the convolution. 

As an interesting by-product, we also find the following fixed-order evolution-equations 
of the first inverse moment of the DAs: 

f\utfhj^= rrf„M5ii«/i+„.l„4(3 + 21n«)l, 
Jo u Jo u [ ni J 

f\ jAu,,l) ^ r.i.M^/l + 2a,ln4f2 + i^')l. (38) 

Jo u Jo u [ l-uj) 

^We use dimensional regularisation with D = 4: + 2e. 
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twist-2: 
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twist-3: 











u(t) 
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V^A (o-part. DAsj 


e 







T (3-part. DA) 


52 




twist-4: 
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chiral-odd 3-part. DAs 


52 





Table 2: Contributions included in the calculation of the correlation functions (26) and 
(27). (5 oc my is the effective expansion parameter of the light-cone expansion; we include 
contributions up to second order in 5 and ccg; those marked by (*) are new. 

These equations allow one to calculate the change of that inverse moment directly for a given 
DA without having to calculate the Gcgcnbauer-momcnts in an intermediate step. The first 
of these relations can also be obtained from the known one-loop evolution kernel of the tt 
twist-2 DA [7], whose anomalous dimensions coincide with those of the relation for 0^ 
is new. 

As we shall see in the next section, the LCSRs do actually not involve the full correlation 
functions, but only their imaginary parts in p^. As in Ref. [1] we take the imaginary part 
only after calculating the convolution integral, which results in closed and comparatively 
simple, albeit lengthy expressions. The distribution amplitudes 0±,||, ^'^^ given by 

their respective conformal expansions, which we truncate at ag. As discussed in Sec. 2, the 
effective expansion parameter of the light-cone expansion is 5, so that the correlation function 
is expanded in both b and ag. We combine both expansions and include terms up to second 
order, i.e. 0(5°'^'^q;°) and 0(5°'^q;^), but drop 0(5^0;^).^° A list of the included terms is given 
in Tab. 2. Note that we have not calculated the radiative corrections to the contributions 
from the 3-particle twist-3 DAs V, A as they are expected to be very small. This follows in 
part from the observation that 0{q.1) terms in the corresponding twist-3 matrix elements do 
also show up in the 0{q1) corrections to g"^""^ and are very small numerically. 

Depending on the specific weak vertex and projection onto the DAs, some diagrams 
contain traces with an odd number of 75, which leads to ambiguities when naive dimen- 
sional regularisation with anticommuting 75 is used. We solve this problem by using Larin's 
prescription for dealing with 75 [18] and replace, whenever necessary, (a^ = CFQ;s/(47r)) 

7/.75 ^ (1 - 4as) ^ e^y^y^y^ I'^^l'^^l'^^ 

75 ^ (1 - 8a,) ^ ey,y,y,y, ^'Y'^'Y^ 

(^tiulb ^ - (1 - OaJ - e^^c./3(7"^. 
^°Terms of 0{5^a1) are not included, either. 
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Note that wc use the Bjorken/Drell convention for the e tensor with 69123 = +1- For the 
special case of the axial- vector form factors and the projection onto the DA g^^\ one can 
implement Larin's prescription by rewriting either the weak vertex or the B vertex. We have 
checked that we obtain the same result in both cases. One might also think of "Larinizing" 
the projection operator onto the DA; the corresponding finite renormalisation will be in- 
dependent due to the nonlocality of the current and is yet unknown. 



4 Numerics 

This section is the heartpiece of our paper, in which we derive the sum rules for S — > y 
form factors and obtain numerical results. The section is organised as follows: in Sec. 4.1 we 
derive the LCSR for one of the seven form factors, V. In Sec. 4.2 we give values for most of 
the needed hadronic input parameters and explain how to determine the sum rule specific 
parameters, i.e. the Borel parameter and the continuum threshold Sq. We also calculate 
Jb^ and Jb^, which are necessary ingredients in the LCSRs. In Sec. 4.3 we motivate the 
need for and introduce models of the twist-2 DAs (j)±^\\. In Sec. 4.4 we calculate the form 
factors at = and discuss their uncertainties. In Sec. 4.5 we present the form factors for 
central input values of the parameters and provide a simple parametrization valid in the full 
kinematical regime of q^. The results for = are collected in Tab. 7, central results for 
arbitrary in Tab. 8. 



4.1 The Sum Rules 

With explicit expressions for the correlation functions in hand, we are now in a position to 
derive the LCSRs for the form factors. Let us choose V{q'^) for a Bq transition as example. 
The corresponding correlation function is Fy as defined in Eq. (26). The basic idea is to 
express Ty in two different ways, as dispersion relation of the expression obtained in light- 
cone expansion on one hand, and as dispersion relation in hadronic contributions on the 
other hand. Equating both representations one obtains a light-cone sum rule for V. One 
side of the equation is hence the light-cone expansion result 



ri>^(pLg')- / ds--^^^^, (39) 

Jml "5 Pb 

with 7cp\P{s,q'^) = Im[ry^], which has to be compared to the physical correlation function 
that also features a cut in p^, starting at m^: 

r'v'%pl,q')= / dsP^^-^; (40) 
the spectral density is given by hadronic contributions and reads 

phys/ 2\ £ 2 2V^(g' ) (-■/ 2 \ I higher— mass states/ 2\ 

Pv {s.q)^fB^mB (5(s-mB)+p+^ (s,g). (41) 

mB + TUv 

Here /^^ is the Bq meson decay constant defined as 

(0|g7M75&|5) = or (m6 + mg)(0|gi756|S) =m|/B,. (42) 
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V 


p UJ K* 6 

r r 


/^[MeV] 

/^(lGcV)[McV] 

/f(2.2GeV)[MeV] 


205 ±9 195 ±3 217 ± 5 231 ±4 
160 ± 10 145 ± 10 170 ± 10 200 ± 10 
147 ±10 133 ±10 156 ±10 183 ± 10 



Table 3: Values of the vector meson couplings, fv is extracted from experiment, fy from 
QCD sum rules for f^/fv, cf. Ref. [19]. 



To obtain a light-cone sum rule for V, one equates the two expressions for Ty and uses 
quark-hadron duality to approximate 

^higher-mass states^^^ ^2) ^ ^LC(^^ ^2^q(^ _ ^^^^ (43) 

where Sq, the so-called continuum threshold is a parameter to be determined within the sum 
rule approach itself. In principle one could now write a sum rule 

and extract V. However, in order to suppress the impact of the approximation (43), one 
subjects both sides of the equation to a Borel transformation 



s — p% ^ ~ Pb 

which ensures that contributions from higher-mass states be sufficiently suppressed and 
improves the convergence of the OPE. We then obtain 



e 



-"^^/^ = r ds e-^/^'p\P{s, q^). (44) 

niB + my J ^2 



This is the final sum rule for V and explains why, as announced in the previous section, 
only the imaginary part of the correlation function is needed. Expressions for the other form 
factors are obtained analogously. The task is now to find sets of parameters (the Borel 
parameter) and Sq (the continuum threshold) such that the resulting form factor does not 
depend too much on the precise values of these parameters; in addition the continuum con- 
tribution, that is the part of the dispersive integral from Sq to oo, which has been subtracted 
from both sides of (44), should not be too large, say less than 30% of the total dispersive 
integral. 



4.2 Hadronic Input Parameters 

After having derived the LCSRs for the form factors, the next step is to fix the parameters 
on which they depend. These are the decay constants of the Bq and Bg meson, /^^ and /b^ , 
the couplings fy"^ of the vector mesons, introduced in Sec. 2, the meson DAs, the quark 
masses mj, and m^, ag and the factorisation scale //ir, and, finally, the sum-rule specific 
parameters and Sq. 
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The fv are known from experiment and are collected in Tab. 3. The fy, on the other 
hand, are not that easily accessible in experiment and hence have to be determined from 
theory. For internal consistency, we determine these parameters from QCD sum rules for 
the ratio fy/ fv, as explained in Ref. [19]. The results are collected in Tab. 3, too. /J had 
already been determined earlier in Ref. [20]; the result agrees with that in Tab. 3. The ratios 
fy/fv have also been determined from lattice [21] and agree with ours within errors. Meson 
DAs are discussed in the next subsection. 

The b quark mass entering our formulas is the one-loop pole mass for which we use 
TOft = (4.80 ± 0.05) GcV, cf. Tab. 6 in Ref. [3]. m^, on the other hand, is the MS running 
mass, ms(2GeV) = 100 Me V, which is an average of two recent lattice determinations [22]; 
the uncertainty in mg has only a minor impact on our results. As for the strong coupling, 
we take as{mz) — 0.118 and use NLO evolution to evaluate it at lower scales. All scale- 
dependent quantities are evaluated at the factorisation scale //ir which separates long- from 
short-distance physics. The only exception arc the form factors Tj, which also depend on 
an ultraviolet scale ^uv which is set to mj. We choose //jr = a/ m?^ ~ ^1 ~ ^-^ GeV as 
reference scale; a variation of /^ir by ±1 GeV has only small impact on the final results. 

The remaining parameters are /b^^, and Sq. fs^^^ has been determined from both 
lattice and QCD sum rule calculations. The state of the art of the former are unquenched 
NRQCD simulations with 2 + 1 hght flavors, yielding fs, = (260 ± 30) McV [23], which is 
slightly larger than the 2003 recommendation fs^ = (240 ± 35) MeV [24]. For /b^, it is 
difficult to find any recent numbers, the consensus being that more calculations at smaller 
quark masses are needed in order to bring the extrapolation to physical rriu^d under sufficient 
control [23]. As for QCD sum rules, both /^^ and /^^ have been determined to 0{af) 
accuracy: fe^ = (208 ± 20) MeV and fs, = (224 ± 21) MeV [25], in agreement with lattice 
determinations. The impact of 0{a^) corrections on /^^^ is nonnegligiblc. As the diagrams 
responsible for these corrections, for instance B vertex corrections, are precisely the same 
that will enter LCSRs at 0{al), we proceed from the assumption that these corrections will 
tend to cancel in the ratio (correlation function) //b. We hence evaluate fs^^ from a QCD 
sum rule to 0{as) accuracy, which reads [26]:^^ 

/so pSQ 



Here {qq) and {qcrgGq) are the the quark and mixed condensate, respectively, for which we 
use the following numerical values at /i — 1 GeV: 

{qq) = -(0.24 ± 0.01)^ GeV^ and (qagGq) = 0.8 GeV^{qq). (46) 

The C are perturbative Wilson coefficients multiplying the condensates. Cqq is known to 
0{as) accuracy [27, 25], CgCq at tree-level. 

The criteria for choosing and Sq in the above sum rule are very similar to those to 
be used for the LCSRs. Ideally, if the correlation function were known exactly, the sum 
rule would be independent of M^. In practice it isn't, but "good" sum rules, plotted as 
function of M^, still exhibit a flat extremum. We hence require the existence of such an 
extremum in and evaluate the sum rule precisely at that point. This eliminates 

^^The contribution of the gluon condensate is not sizable and we therefore neglect it. 
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So fe 


So M2 fB 


sctl 
set 2 
sets 


4.85 
4.80 
4.75 


33.8 3.8 148 
34.2 4.1 161 
34.6 4.4 174 


34.9 4.2 169 
35.4 4.4 183 
35.9 4.6 197 



Table 4: Parameter sets for fs^ and Jb^ to 0{as) accuracy, fs^ and fs^ are given in MeV, 
So and in GeV^. Note that the values of fs^ , given in the table are not to be interpreted 
as meaningful determinations of these quantities, cf. text. 

as independent parameter and leaves us with sq. As already mentioned after Eq. (44), the 
purpose of the Borel transformation is to enhance the contribution of the ground state to 
the physical spectral function with respect to that of higher states. We hence require that 
that continuum contribution, that is the integral over for s > sq, must not be too 

large. To be specific, we require 

(^j^ dsp"'\s)e-''^'^ I (^j^ dsp''''{s)e-'/'^'^ < 30%. 

This puts a lower bound on sq. The larger sq, the smaller M^, the position of the minimum, 
and the larger nonperturbative contributions to (45). As the condensates are meant to yield 
small nonperturbative corrections, but blow up at small M^, requiring the nonperturbative 
corrections to be not too large puts an upper bound on sq. For fs^^, we require the highest 
term in the condensate expansion, the mixed condensate, to contribute less than 10% to the 
correlation function. For LCSRs, which rely on an expansion in higher twist rather than 
higher condensates, we correspondingly require the contribution of higher twists to the LCSR 
not to exceed 10%. One more requirement on the Sq is that they not stray away too much 
from "reasonable" values: sq is to separate the ground state from higher mass contributions, 
and hence should be below the next known clear resonance in that channel. Assuming an 
excitation energy of 0.4 to 0.8 GeV, we thus expect the sq to lie in the interval 32 to 37 GeV^, 
which is evidently fulfilled by all sq quoted in Tab. 5. 

Applying the above criteria to (45), we obtain the sets of (so,M^) collected in Tab. 4, 
together with the resulting /b^ ^,. We would like to stress that these values are not to be 
interpreted as new independent determinations of fsq^^j but are intermediate results to be 
used in the evaluation of the LCSRs. 

We proceed to determine the continuum thresholds and Borel parameters for the LCSRs, 
using the same criteria as above. In order to keep the complexity of the calculation at a 
manageable level, for each form factor the corresponding set is determined only once, at 

— 0. To avoid confusion between parameters entering (45) and those entering the LCSRs, 
let us call the latter ones MIq and Sq where F is the form factor. For larger q'^, these 
parameters are expected to change slightly. Part of this effect can be taken into account in 
the following way: the tree-level LCSR to twist-2 accuracy reads, basically. 
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rrib 


4 

U C 


U C 


4' cf' 

U C 


4' 4' 

U C 


4 

U C 


So' 

u c 


set 1 
set 2 
set 3 


4.85 
4.80 
4.75 


35.2 1.7 
35.8 2.1 
36.4 2.1 


33.0 1.7 
33.6 1.6 
34.2 1.6 


33.7 1.7 
34.2 1.8 
34.7 1.9 


34.1 1.7 
34.7 1.8 
35.3 1.9 


34.8 1.7 
35.3 1.9 
35.8 2.1 


34.7 1.7 
35.2 1.8 
35.7 1.9 



Table 5: Parameter sets for Bg ^ p for V, Aq, Ai, A2, Ti and T3. As Ti(0) = T2(0) the 
corresponding parameters are equal. Sq and in GeV^. 



which implies that the expansion parameter is uMIq rather than M^q. We hence rescale 
the Borel parameter as 

Mlc^MlJ{u){q^) 
with the average value of rt, {u){(f), given by 



\Juo U 



-{ml-il-u)q^)/{uMlc) 



2\ n Qd „„A lr.,\l^ A r~'„'\r2 



with, approximately, (m)(OGcV ) = 0.86 and (■u)(14GeV ) = 0.77. The optimum Borel pa- 
rameter hence becomes larger with increasing g^, which agrees with what one finds when M^^q 
is determined without rescaling. Parametrising the relation between the Borel parameters 
of local and hght-cone correlation functions as 



we obtain, for 
similar. 



M^c = CeMV(«), (47) 
p, the values collected in Tab. 5. The sets for other transitions are 



4.3 Models for Distribution Amplitudes 

As mentioned in Sec. 2 and detailed in App. B, the DAs entering the LCSRs can be modelled 
by a truncated conformal expansion. It turns out that the dominant contributions to the sum 
rules come from the twist-2 DAs 0±,||, which to NLO in the conformal expansion are described 
by the lowest three Gegenbauer moments: aj" '" = 1, which follows from the normalisation 
of the DAs, af'", which is nonzero only for i^*, and a^'". In Ref. [4], it was these three 
parameters that were used to define the models for all terms a„>3 were dropped. 

The numerical values of a\^i (and higher moments) are largely unknown. a\ has been 
determined from QGD sum rules in [28, 19, 29]. Averaging over the determinations, we 
choose 

o}\{K\ 1 GeV) = 0.10 ± 0.07 = a\{K\ 1 GeV) (48) 

as our preferred values. Note that positive a\ refer to a K* containing an s quark - for a K* 
with an s quark, ax changes sign. 

Predictions for 02 '" also come from QGD sum rules [14, 28, 19, 20] and read 

a|(p, 1 GeV) = 0.18 ±0.10, a5(i^*, 1 GeV) = 0.09 ± 0.05, a|(0, 1 GeV) = ± 0.1, 

a^(p,lGeV) = 0.2±0.1, a^(X*, 1 GeV) = 0.13 ± 0.08, 0^(0, 1 GeV) = ± 0.1. 

(49) 
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All these determinations have to be taken cum grano salis, as the sum rules do not exhibit 
a clear Borel-window and also become increasingly unreliable for larger n}^ 

But even assuming ai,2 were known to sufficient accuracy - under what conditions is a 
truncation of d after 02 is justified? We have seen in Sec. 2 that after the convolution with 
a smooth short-distance function T the contributions of higher a„ fall off sharply. So the 
actual question is not so much how the truncated expansion compares to the full convolution 
integral, but rather how the neglected terms compare to other terms, for instance originating 
from 3-particle DAs, which are included in the LCSR. For instance, assuming Oj > 0.05 it 
is necessary to include a| and 02^4^6,8 order to match the size of the contributions from 
quark-quark-gluon matrix elements, and even for aj > 0.01 one still needs a| and 0^4. If, 
on the other hand, one consistently neglects terms that contribute less than 1% to the form 
factor, one can drop nearly all contributions from quark-quark-gluon matrix elements, unless 
their values as given in App. B are grossly underestimated. If > 0.05, one then has to keep 
a^4 g, but can drop all ai[>o- The upshot is that, in view of the lack of information on On , 
it is a good idea to devise models for 0^ y with a small number of parameters, possibly tied 
to experimental obscrvables, and a well-defined "tail" of higher-order Gegenbauer moments. 
This task is undertaken in Ref. [17]. 

Following Ref. [17], we introduce two-parameter models for (the symmetric part of) (f)±^\\ 
which are defined by the fall-off behavior of the Gegenbauer moments a„ in n and the value 
of the integral 

du^ = l + J2(^2n; (50) 

n=l 

A = 1 for the asymptotic DA. In particular we require A to be finite, which implies that 
the a„ must fall off sufficiently fast. The choice of A as characteristic parameter of (f) relies 
on the fact that it is directly related to an experimental observable, at least for tt and rj, 
namely the 7r(?7)77* transition form factor, for which experimental constraints exist from 
CLEO [30] . We assume that the vector meson DAs are not fundamentally different and take 
the range of A(7r) extracted from CLEO as the likely range for A(p). The second parameter 
characterizing our models is the fall-off behavior of the Gegenbauer moments in n, which 
we assume to be powerlike. We then can define a model DA 0+ in terms of its Gegenbauer 
moments 

^ (51) 



" (n/2 + l)« • 
using the generating function of the Gegenbauer-polynomials, 



00 



this model can be summed to all orders in the Gegenbauer-expansion: 

4>Uu) = I' dt{- Inty-' (/(2« f{2u - 1, -^/^)) . (52) 
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This is due to the different power behavior of perturbative and nonperturbative terms in n, cf. Ref. [20] . 
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Figure 2: Examples for model DAs 0+ as functions of u, for A = 1.2 and a = 1.5,2,3,4 
(solid curves), as compared to the asymptotic DA (dashed curve). For a ^ 1, 4>^ approaches 
the asymptotic DA. 

0.2 
0.1 
0. 
-0.1 

0. 0.2 0.4 0.6 0.8 1. 

Figure 3: Models for the asymmetric contributions to the twist-2 DA for ai = 0.1. Sohd 
curves: 0^^™'"*" as function of w for 6 e {2, 3, 5}; dashed curve: 0^^™'+. 

The corresponding value of A is A+ = ({a). In order to obtain models for arbitrary values 
of A, we split off the asymptotic DA and write 

0+(A) = 6uu + _ Quu) . (53) 

Evidently one recovers the asymptotic DA for A = 1 and the truncated conformal expansion 
with 02 = A — 1 and an>4 = for a ^ oo. The above formula is only valid for a > 1, as 
otherwise A+ diverges, or, equivalently, 0^ does not vanish at the endpoints u = 0,1. In 
Fig. 2 we plot several examples of 0^ for a fixed value of A. 

Our preferred values for A, a and the corresponding values of (1 GeV) are collected 
in Tab. 6. We choose a = 3 ± 1 in order to obtain nonneghgible effects from higher order a„. 
The choice of A(p) is motivated by the fact that all available calculations indicate 02 > 0, 
hence A > 1. We then fix the maximum A in such a way that it yields 02 < 0.2, which, 
given the fact that the sum rule results (49) arc likely to overshoot the true value of 02, 
appears as to be the likely maximum value. We then obtain A(p) = 1.15 ± 0.10, with a 
rather conservative error. We choose the same values for cv. For K* and 0, we take into 
account that the values of 02 appear to have the tendency to decrease, which was noticed 
already in Ref. [28]. Assuming that the decrease is 20% from p to K*, and another 20% 
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A a 


„-L,|l/'-| r^T^T\ „-'-i! /I 

02 [i- (jreV ) [i LreV ) 


K* 


1.15 ±0.10 3±1 
1.12 ±0.10 3±1 
1.10 ±0.10 3±1 


Omtofj 0.03 ±0.02 
0.07l[J:[J? 0.02 ±0.01 
06+°-°^ 02+° °^ 



Table 6: Values for A(l GeV) = A-*-'" (1 GeV), a and the corresponding values of 03^4 (1 GeV). 

from K* to (f), we arrive at the numbers quoted in Tab. 6. 

Evidently, DAs defined in dependence of A also require a specification of the scale at 
which they are valid. As we presume that A(//) will be measured, if at all, at the low scale 
/i = 1 GeV, we choose this as the reference scale. A at higher scales can be obtained from 
Eq. (50), using the leading-order RG-improved expressions for a„(//), or, if /j, is not too 
different from 1 GeV, from the unimproved expression Eq. (38). 

Models for the asymmetric part of the DA, relevant for K*, can be constructed in a 
similar way as 

i^tH = ffy ^^(- 1"^)'"' (/(2^^ - 1, v^) - mu - 1, -v^)) . (54) 

One relevant parameter is b, and as the second one we choose Oi. Models for the asymmetric 
part of 4> with arbitrary oi can then be defined as 

C^'^'+ = ai(3/2)^Vi,+ (x.). (55) 

Examples for such models are shown in Fig. 3. 
0asym,+ ^-^^ contributes to the value of A: 

^a.ym,+ ^ / l"^ ^ -Oi (3/2)^C(^ 3/2), 

Jo 

where C(^) — YlT=o^/(^ + Hurwitz ( function. Our models for the K* DA 

are hence characterized by four parameters: A, a of the symmetric part and Oi, b of the 
asymmetric part. The total value of A is given by 

^total,+ — ^ _|_ ^asym,+ 

In the actual calculation we choose a — b. 
4.4 Results for = 

Let us first analyse the form factors for = 0. Using the input parameters given in Tabs. 5 
and 6, we obtain the results collected in Tab. 7. 

For the discussion of theoretical uncertainties, we distinguish between uncertainties that 
can be reduced by future more accurate determinations of the corresponding hadronic pa- 
rameters and others that are either systematic uncertainties, inherent to the method of 
LCSRs, or parameter uncertainties not likely to be reduced in the near future. The latter 



22 





m 


Am, 






At 


Atot 


A„, 






0.323 


0.007 


0.025 


0.005 


0.013 


0.029 




. R - 
^0 




0.303 


0.004 


0.026 


0.009 


0.006 


0.028 




- R 

Af' 




0.242 


0.007 


0.020 


0.004 


0.010 


0.024 




Ap 


~*p 


0.221 


0.008 


0.018 


0.002 


0.011 


0.023 








0.267 


0.004 


0.018 


0.004 


0.010 


0.021 




IJlBq- 

^3 




0.176 


0.001 


0.013 


0.001 


0.009 


0.016 




yBs- 




0.311 


0.006 


0.021 


0.003 


0.013 


0.026 


-0.43 5a, 


aB,- 
^0 




0.363 


0.003 


0.032 


0.006 


0.009 


0.034 


-0.37 


Af^- 


*K* 


0.233 


0.007 


0.019 


0.002 


0.010 


0.023 


-0.32(5,, 


aBs-^k* 


0.181 


0.008 


0.021 


0.001 


0.010 


0.025 


-0.30 5,,, 


rpBs-^K* 


0.260 


0.005 


0.021 


0.003 


0.010 


0.024 


-0.33(5,, 


rpBs^K* 

^3 


0.136 


0.003 


0.013 


0.000 


0.008 


0.016 


-0.175,, 




*K* 


0.411 


0.008 


0.029 


0.003 


0.013 


0.033 


0.44 5,, 


A^q- 
^0 


*K* 


0.374 


0.009 


0.031 


0.005 


0.008 


0.034 


0.39 5„, 


Af"- 


*K* 


0.292 


0.009 


0.025 


0.002 


0.009 


0.028 


0.33 5a, 




*K* 


0.259 


0.009 


0.023 


0.001 


0.010 


0.027 


0.31 5„, 




*K* 


0.333 


0.005 


0.026 


0.003 


0.010 


0.028 


0.34 5„, 


rpBq — 

^3 


^K* 


0.202 


0.002 


0.016 


0.001 


0.008 


0.018 


0.18 5,, 




->a; 


0.293 


0.006 


0.025 


0.002 


0.013 


0.029 




^0 


-^LU 


0.281 


0.012 


0.027 


0.003 


0.006 


0.030 




Af'- 




0.219 


0.008 


0.021 


0.001 


0.010 


0.025 






-^LU 


0.198 


0.007 


0.018 


0.001 


0.011 


0.022 




rj^Bq- 


-^LU 


0.242 


0.003 


0.019 


0.002 


0.010 


0.022 




rpBq- 

^3 


-^LJ 


0.155 


0.000 


0.012 


0.000 


0.009 


0.015 




yBs- 


-^4> 


0.434 


0.004 


0.032 


0.003 


0.014 


0.035 




Af- 


-0 


0.474 


0.002 


0.031 


0.005 


0.019 


0.037 






-0 


0.311 


0.007 


0.027 


0.002 


0.009 


0.029 




A^ 


-^4> 


0.234 


0.011 


0.024 


0.001 


0.009 


0.028 




rj^Bs- 




0.349 


0.004 


0.031 


0.002 


0.010 


0.033 




rpBs- 




0.175 


0.003 


0.016 


0.000 


0.007 


0.018 





Table 7: Form factors at = for parameter set 2 of Tabs. 4 and 5, i.e. = 4.8 GeV. The 
form factors are defined in Eqs. (2) and (4). The penguin form factors Tj are evaluated at 
the UV scale = mi,. is the variation of the result with mt, i.e. the maximum deviation 
between the results obtained for sets 1, 2 and 3. Ayp is the maximum deviation found by 
scanning the 7-parameter space discussed in the text. A^ and A^ are the uncertainties 
induced by the vector and tensor couphngs in Tab. 3. The total error Atot is obtained by 
adding A(mi,,7p,L,T) in quadrature. Form factors involving K* carry one more uncertainty 
A„ induced by the Gegenbauer moment oi, with 5„ = [ai{K* , 1 GeV) — 0.1]. 
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comprise the dependence of the form factors on the LCSR parameters sq, M^, and, via 
Jb, the quark and mixed condensate. Our results also depend, very mildly, on and, 
more importantly, on the meson DAs which are described by the 2-parameter model (53). 
All these parameters induce a theoretical error of the form factors which we determine by 
varying 

• the threshold Sq by ±1.0 GeV^; 

• the Borel parameter in Eq. (47) by ±1.5 GeV^; 

• the infrared factorization scale //jr = •>/ — ml by ±1 GeV; 

• the quark condensate and the mixed condensate as indicated in Eq. (46); 

• the first inverse moment of the twist-2 DAs, A, by ±0.1; 

• the power behavior of the Gegenbauer moments, a, by ±1; 

• the strange quark mass by ±20%. 

The largest deviation of the form factor from its central value, in this 7-parameter space, is 
dubbed Ay^ and amounts to typically 7 to 11%. In Fig. 4 we show the dependence of selected 
form factors on A and a. The uncertainty in these parameters is the most important single 
source of error of the form factors and amounts to half of the total error. 

The form factors also depend, rather mildly, on mb'. varying by ±0.05 GcV around 
the central value 4.8 GeV, and using sq and as given in Tabs. 4 and 5, we obtain the 
error A^^ which ranges from 1% to 5%. 

One more source of uncertainty of the form factors is due to fv and fy, the vector and 
tensor coupling of the vector mesons. This is easily understood by splitting the generic form 
factor F into two terms proportional to and = f: 

F = f^F^ + fF'^. (56) 

As argued in Sec. 2, the first term is of order 5 oc my^ the second of order 1 and indeed, 
for most form factors, is the dominant contribution. The present errors of as collected 
in Tab. 3, are nonneghgible. is accessible to lattice calculations and first results have 
been reported in Ref. [21], which indicates that a reduction of the error of seems feasible. 
In order to allow the adjustment of our form factors to new results for /-^, we give explicit 
results for F^ and F^ in App. A. The uncertainties ^t,l of the form factors due to the 
present values of /^'^ are included in Tab. 7. At is typically of order 4%, A^ is much 
smaller. 

For transitions involving the K* , an additional uncertainty is induced by the first Gegen- 
bauer moment ai, and is given by A^-^ in Tab. 7, where the quantity is defined as 
[ai(i^r*, 1 GeV) — 0.1]. Note that ai{K*) refers to a sq bound state and hence ai{K*) = 
—ai{K*), which explains the negative sign of the corresponding entries in Tab. 7. Again we 
aim to make our results adjustable to any future improvement in the determination of oi 
and give explicit results for the corresponding contributions in App. A. 

Some important features of the results collected in Tab. 7 are: 
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Figure 4: The form factors ^(0), y4i(0), ^2(0) for Bq ^ p and ri(0) for ^ K* as 
functions of A, the main parameter of the twist-2 DAs. Sohd hues: central values of input 
parameters. Dashed lines: variation of the form factors with a change of a, the second 
parameter of the DAs, by ±1. Allowed values of A: cf. Tab. 6. 



• the form factors for Bq — > K* transition are about 20% larger than those for Bq — > p. 
The reason for this is twofold: on the one hand, the K* vector and tensor couphngs 
are larger than those of the p. On the other hand, the SU(3)-breaking of the twist-2 
DAs, parametrised by the first Gegenbauer moment ai, gives a positive contribution 
to the form factors; 

• the form factors for Bs — > K* have a tendency to be smaller than those for Bq p. 
The reason for this is a negative contribution of ai and the fact that /^^ is larger 
than fBq- On the other hand, the optimum Sq are also larger than for Bq — > p, which 
partially compensates the first two effects; 

• the Bq ^ u form factors are slightly smaller than those for Bq p. This is a 
consequence of the fact that the uj vector and tensor couplings are smaller than those 
of the p; 

• the total theoretical error is dominated by that of the twist-2 DAs and the sum rule 
parameters sq and M^. The former can, in principle, be reduced by future calculations, 
the second is systematic and irreducible. 

The typical total uncertainty of each form factor is 10%, ranging between 8% and 13%. 
Any significant reduction of the error requires more accurate information on the twist-2 DAs. 
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The minimum irreducible theoretical uncertainty is set by the systematic uncertainty of the 
LCSR approach and encoded in the dependence of the results on sq and M^; it amounts to 
about 6 to 7%. 

Let us also compare our results to those obtained in Ref. [4] by the same method, but 
with less sophistication. The main difference between our present and our previous analysis 
is the inclusion of radiative corrections to 2-particle twist-3 contributions, the description of 
the twist-2 DAs by models including all-order effects in the conformal expansion, the more 
accurate determination of the sum rule parameters Sq and and the much more detailed 
error analysis. The most striking difference between the actual results affects the Bg K* 
transition, whose form factors were predicted, in [4], to be between 10 and 30% smaller than 
those in Tab. 7. The reason for this discrepancy is mainly the more accurate determination 
of and sq we employ in the present analysis — in Ref. [4] all form factors were determined 
for the same values of and sq- All other form factors quoted in [4] agree, within ±15%, 
with those of Tab. 7, which is within the theoretical uncertainty stated in [4]. The only 
exception is T^,{^), which deviates by between 15 and 45% from the numbers obtained in 
[4]. The reason for this discrepancy lies in the (correct) treatment of factors \/{pz) in our 
present paper, cf. Sec. 3. 

4.5 Results for 7^ 0, Fits and Extrapolations 

In this subsection we discuss the g^-dependence of the form factors. The results of the 
LCSR calculation are plotted in Figs. 5 and 6. They can be parametrised in terms of simple 
formulas with 2 or 3 parameters, which are valid in the full kinematical regime in (f. The 
corresponding parameters are collected in Tab. 8. 

As mentioned in Sec. 2, the LCSR method is valid for large energies of the final state 
vector meson, Ey ^ Aqcd, which implies, via the relation = m\ — 2mBEv, a restriction 
to not too large q'^. We include values in the regime 



which has to be compared with the maximum physical 9phys,max — {^b — rnvY of 20.3 GeV^ 
for Bq (p, u), 19.2 GeV^ for B^^K*, 20.0 GeV^ for ^ K* and 18.2 GeV^ for ^ 0. 
The main aim of this subsection is to provide fits for the LCSR results, which are valid in 
the full physical regime of q^. We will comment below on the dependence of the fit results 
on the actual value used for ^lcsr max- 

We closely follow the procedure we used in our previous paper on S — > pseudoscalar form 
factors, Ref. [1] . Generically, barring the occurence of anomalous thresholds, any form factor 
F(g^) has singularities (poles and cuts) for positive real g^, starting at the position of the 
lightest resonance coupling to the relevant current, and hence can be written as a dispersion 
integral in q^. Splitting off the lowest-lying resonance with mass ruR, one has 



where is the threshold for multiparticle contributions, which can be above or below m\. 
Keeping only the first term and neglecting the integral altogether one obtains the vector 



0<g^<gEcsR,max = 14GeV% 



(57) 




(58) 
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Figure 5: Form factors for Bg decays as functions of q^, for central values of input parameters. 

meson dominance (VMD) approximation.^^ Even though this approximation is expected to 
work very well close to the pole, it certainly won't work far away from it, e.g. at = 0. For 
S — > TT transitions it was argued in Ref. [31] that the integral can be modelled by a second 
pole at larger q^, which is unrelated to any physical resonance: 

with the three independent parameters ri^2 and mfit. 

The dominant poles at q^ — m\ correspond to resonances with quantum numbers — 
1' for V and Ti, 0~ for and 1+ for ^1^2,3 and T2^3, 2^. As discussed in Sec. 2, not all these 

^^This notion conies from the analysis of electromagnetic form factors, where the first resonance is the p. 
In weak decays, however, the lowest resonance is, in general, not a vector meson, so that the notion VMD 
is, strictly speaking, obsolete. 
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Table 8: Fits for the form factors valid for general q^. Columns 2 to 4 give the results of 
Tab. 7 for = 0, including the errors Atot and A^, . The remaining colums give the fit 
parameters. Note that we fit the form factor T3, defined in Eq. (10), instead of T3. The 
fit formulas to use are given in the last column, the masses uir are given in Tab. 9. The 
penultimate column gives the fit error 5 as defined in Eq. (62). 
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Figure 6: Form factors for Bg decays as functions of g^, for central values of input parameters. 
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5.28 5.63 5.32 5.68 
5.37 5.72 5.42 5.77 



Table 9: B meson masses in units GeV, taken from Ref. [32]. 



form factors are independent, and the question arises which ones to fit to the above equation 
- or any similar formula - and which ones to define in terms of the others. As Eq. (59) 
contains two explicit poles, we decide the above question in favor of the form factors with 
the steepest increase in q^, which means that the independent form factors are V, ^0,1,2 and 
7i,2, Ts, whereas T3 and A3 are the dependent ones, defined as in Eqs. (3) and (8). 

The values of the resonance masses mn in (59) are known from experiment for 0~ and 
1~ in the Bg channel and 0~ in the Bg channel; the other masses are obtained using heavy 
quark symmetry relations [32], the numerical values are collected in Tab. 9. 

We shall use fits to Eq. (59) for the form factors V, Aq and Ti, where the lowest pole 
m|j lies well below the multiparticle threshold {mBg^^ + ttIt^^kY. If, on the other hand, the 
lowest physical pole lies sufficiently close to the multiparticle threshold or even above it, 
then it may be impossible to "resolve" the poles from a low-g^ "perspective". In this case it 
is more appropriate to expand the form factor to second order around the pole, yielding 



Fie) - 



ri 



+ 



r2 



fit 



(60) 



with the three parameters ri 2 and mfif This is the fit formula we shall use for the axialvector 
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form factors, in particular A2 and T3. For Ai and T2, on the other hand, the residue of the 
double pole in mat turns out to be extremely small, so that it can be dropped and one is 
back to the VMD formula 

Hq') = 1 2 1 (61) 

albeit with an effective pole mass met unrelated to any resonance. 

The fits of the LCSR results to the above formulas are collected in Tab. 8; they differ 
from the LCSR results obtained for < 14GeV^, by no more than 2.5%. In Tab. 8 we 
indicate the "quality" of the fit by 6, which is the maximum deviation of the fit relative to 
the mean value of the form factor in percent and defined as 

where the sum runs over t e {0, 0.5, 1, . . . , 14}. 

We have also tried fits to the two pole ansatz (59) without fixing one of the masses. In 
this case the lowest pole is fitted to lie below the actual resonance pole, by up to 1.5 GeV^. 
Given the fact that LCSRs are valid for small far away from the pole, one cannot expect 
them to resolve its position with perfect accuracy. Nonetheless we take it as an indication 
for the consistency of our approach that the double pole formula with unrestricted pole 
positions gives results that agree quahtatively with those from the restricted fits. We also 
have checked the dependence of the fits on the maximum value of ^maxLCSR ^'^ which 
LCSR results are included into the fit. It turns out that the fits are very robust against 
lowering ^maxLCSR' lowering it from 14GeV^ to 7GeV^ changes the fitted values at 20 GeV^ 
by at most 8%, T2 being the odd one sticking out. In Fig. 7 we show the effects of a change 
of ?^ax,LCSR on T.^^" and Af-"". 

Let us now turn to a consistency check of our fits. One can express the residues of V, Ti 
and Aq for B ^ p in terms of decay constants and strong couplings as follows: 



^1=^^ fB*9BB*p, r^' = —gBB*p, r^" = 77— 9BBp, (63) 



^^r_Y_^rns + rrwf^^^^^ (64) 



where is the tensor coupling of the B* meson defined in the same way as light vector 
tensor couplings, Eq. (21). /b has been discussed in Sec. 4.2; its value is about 200 MeV 

and we expect Jb* and /J, to be of about the same size. The values of the strong couplings 
gBBp and Qbb'p are more controversial as discussed below. As a first check, consider the 
gf-independent ratio 

_ '1 _ 

rfi niB* fs* 

The fitted values of ri are collected in Tab. 10 and yield agt = 1.16 — very close to (64). 
For ri fitted using parameter sets 1 and 3 we find a — 1.16 and 1.17, respectively. Any 
further check of the fitted ri requires information on the couplings gBB(*)pj which have been 
calculated from both LCSRs [33] and within the constituent quark meson (CQM) model [34] 
— with significantly different results. The situation resembles that for gDD*-^-, where LCSR 
determinations are typically by a factor 2 smaller than lattice and CQM calculations [35]. For 
this coupling there actually exists an experimental measurement by CLEO [36], which agrees 
with the lattice and CQM determinations, but disagrees with LCSRs. For the corresponding 
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0. 5. 10. 15. 20. 0. 5. 10. 15. 20. 

Figure 7: Comparison of the consistency of fits of Ti ' and '' obtained for different 
values of g'LcsRmax- Dots: LCSR results for < 14GeV^. Lines: fits according to Eqs. (59) 
for Ti and (61) for Ai, for Q'LcsRmax between 7 and 14GeV^. The maximum discrepancy 
between the fit-results at q'^ = 20GeV^ is 2% for Ti and 5% for Ai. 

B coupling Qbe'tt there is no experimental measurement, as the decay B* — > Stt is forbidden 
by phase space, but one can use heavy quark scaling to obtain gBB*-n from the measured 
gDD*TT and compare it with the corresponding theoretical predictions. It turns out that again 
lattice and CQM calculations are favored, whereas the LCSR calculation gives a too small 
result, which can be understood following the discussion in Ref. [37]. The recent LCSR 
determination [33] has up-to-date input parameters and they get from a tree-level analysis 

9BBp = 5.37, gBB'p = 5.70 GeV"^ . (65) 

For pseudoscalar mesons, NLO calculations have consistently yielded smaller values than 
tree level determinations, cf. Rcf. [38], which, if true also for the p. would widen the gap 
between the results from different methods even further. The CQM-model predictions are 
[34]: 

9BBp = -V2/?^ = 7.2 with /? = -0.86, (66) 
gBB*p = VSX . I^a^ ^ 10.0 GeV-^ with A = 0.6 GeV^^ (67) 

V "^B Jvr 

It is hard for us to judge on the validity of this approach, but as far as we understand the 
model is further based on empirical success. In Tab. 10 we compare the residues for the 
S — > p transition as obtained from our fits. Tab. 4, to their values given in Eq. (63), using 
the couphngs (65), (66) and (67). For V and Ti with a 1~ pole the CQM residues are about 
10% larger, and the LCSR about 40% lower than the fitted values. As discussed above, the 
LCSR results are expected to fall short of the real values, so this is an excellent confirmation 
of our results. The Aq form factor shows some discrepancy which may indicate that either the 
estimate of the Qbbp coupling is too low or that the second pole in the fit, ~ 33 GeV^, 
is too close to the resonance pole to allow a clean determination of its residue. Taken 
altogether, however, the agreement of our fitted results to that of independent calculations 
is an excellent confirmation or our results. 
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our nt 


LUoK [ooj 


UL^M [o4J 


'1 


1.05 


0.65 


1.14 


' I 


0.90 


0.57 


1.00 


r-Ai) 
'1 


1.53 


0.70 


0.94 



Table 10: Residues of the lowest-lying pole for V^?"*'', and A^''^'' obtained from our 

fits as compared to Eq. (63) with input values from LCSR and CQM determinations. 

5 Summary & Conclusions 

In this paper we present a thorough and careful examination of the predictions of QCD sum 
rules on the light-cone for the form factors of Bq and Bg transitions to p, a;, K* and 0. Our 
main results for zero momentum transfer = are collected in Tab. 7, those for general 
in Tab. 8. 

The present analysis is a sequel of our work on B pseudoscalar form factors, Ref. [1], 
and an extension of the previous work of one of us on ^ vector form factors, Ref. [4] . It 
improves upon the latter by 

• including predictions for all form factors of Bq g transitions to 0{as) accuracy for twist- 
2 and 3 2-particle contributions; 

• a more sophisticated method for fixing sum rule specific parameters, cf. Sec. 4.2; 

• implementing recently developed new models for the dominant nonperturbative vector 
meson contributions, the twist-2 vector meson distribution amplitudes, cf. Ref. [17]; 

• the possibility to implement future updates of some hadronic parameters in a straight- 
forward way, cf. App. A; 

• a careful assessment of uncertainties at zero momentum transfer, cf. Sec. 4.4; 

• a parametrization of the g^-depcndcnce of form factors valid in the full physical regime 
of momentum transfer, cf. Sec. 4.5; 

• a variety of consistency checks for the robustness of the g'^-fits and their numerical 
results. 

The accuracy of our results is limited, on the one hand, by the uncertainty of hadronic 
input parameters and, on the other hand, by the systematic uncertainty induced by the fact 
that QCD sum rules on the light-cone are an approximative method. The uncertainty due to 
the variation of only the sum rule specific parameters is about 7%, which cannot be reduced 
any further and hence sets the minimum theoretical uncertainty that can be achieved within 
this method. An equally large theoretical uncertainty is induced by hadronic parameters and 
can, in principle, be improved upon. We quote in particular the tensor couplings fy of vector 
mesons, which presently come with the rather large error quoted in Tab. 3. Improvement 
should be possible by dedicated lattice calculations, a first example of which is Ref. [21]. 
Another relevant hadronic parameter is A, the first inverse moment of the twist-2 vector 
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meson distribution amplitudes, as defined in Sec. 4.3. We have inferred a likely range for this 
parameter for p and uj mesons from the known experimental constraints on A(7r), and further 
determined a range for A(ir*) and A(0) from the observed decrease, within QCD sum rule 
calculations, of the second Gcgcnbaucr moment 02 with increasing meson mass. Comparing 
the theoretical errors collected in Tab. 7 with the global theoretical uncertainty ~ 15% quoted 
in our previous publication [4], we have achieved a reduction to about 10%. This is partially 
due to a reduction of the uncertainties of the hadronic input parameters, in particular m;,, 
and partially due to a refinement of the assessment of sum rule specific uncertainties as 
discussed in Sec. 4.2. Any future reduction of the total uncertainty will depend on more 
accurate determinations of A, which are absolutely essential not only for light-cone sum 
rule calculations, but also for exploiting the full potential of QCD factorisation formulas for 
nonleptonic exclusive B decays [6]. We take this occasion to urge lattice practitioners to 
take up the challenge and develop new and ingenious methods to tackle this problem — or 
just give us an accurate value of which would already be a big step forward. 

The prospects for future direct determinations oi B ^ V form factors from lattice cal- 
culations do appear a bit clouded. On the one hand, there are two recent studies, by the 
SPQcdR and UKQCD collaborations, Ref. [40], using an improved Wilson action and the 
quenched approximation. The h quarks are fully relativistic and have typical masses of about 
2 to 3 GeV, so they need to be extrapolated to the physical h quark mass. On the other 
hand, we conclude from [39] that an unquenched calculation in NRQCD is not really on the 
menu, which, as far as we understand, is due to an improvement in the treatment of light 
quark masses on the lattice, causing the p and other vector mesons to become instable par- 
ticles without a pronounced plateau in the fall-off of the correlation function, and essentially 
prevents a precise determination of their properties from lattice. We do not pretend to be 
sufficient experts in LQCD to be able to meaningfully comment on these issues, but remain 
hopeful that the situation will be clarified in due course. 

We have calculated all form factors for < < 14GeV^; the upper bound on is due 
to the limitations of the light-cone expansion which requires the final-state meson to have 
energies E ^ Aqcd: for q^^x = 14GeV^ the meson energy is E = 1.3 GeV. In order to 
facilitate the use of our results we have given, in Sec. 4.5, Eqs. (59), (60) and (61), simple 
parametrisations that include the main features of the analytical properties of the form 
factors and are valid in the full physical regime < < (m^ — myY- The corresponding 
results for our preferred set of input parameters are given in Tab. 8. We have checked that 
the fit results arc fairly insensitive to the maximum value of included, and that reducing 
the latter to e.g. 7 GeV^ changes the extrapolated values of the form factors at q^ — 20 GeV^ 
by typically only 1 to 2%, and by 8% at most (for T2). 

In Sec. 1 we mentioned factorisation formulas for form factors derived in SCET, Ref. [8, 9, 
10], which in particular imply certain (heavy quark) symmetry relations. Since the objective 
of this paper was to provide numerical results, ready for use in phenomenological applications, 
we did not discuss the question whether and to what extent our results fulfill these relations, 
nor the size of symmetry-breaking corrections. A previous study of the corresponding effect 
in 5 — i> pseudoscalar decays has indicated that such corrections are likely to be nonnegligible 
[41]. We plan to come back to these points in a future publication. 
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Appendix 

A Form Factors For Different fy^ and ai{K*) 

The form factors can be written as sum of two contributions which arc proportional to the 
vector meson's vector coupling fv = and the tensor coupling fy = f'^, respectively. The 
uncertainties of these parameters as tabled in Tab. 3 are nonnegligible, but amenable to fu- 
ture improvement by e.g. lattice calculations, cf. Ref. [21]. The same applies to ai{K*) which 
also comes with a considerable uncertainty, cf. Eq. (48). In order to allow the adjustment of 
our results to improved determinations of these parameters, we write the generic form factor 
F as 

F = /^(F^ + af + f^{F^ + a[ , (A.l) 

where the hatted quantities are normalized to the central values used in our calculation, 
i.e. the couplings of Tab. 3 and af = a^, = a| as given in (48). For instance fp — 
/p/(205MeV). Note that ai{K*) = ai{K*) and that hatted quantities are trivially invariant 
under LO scaling. and are collected, for = 0, in Tab. A and F"^ in Tab. B. To 
give an example, for Aq"~*^ we obtain 

A^'-^^UO) = ( 0.2469- ^OL^ 0.0398^ 

° ^ ^ 217MeV V 0.1 J 

flJl GeV) / af*(l GeV) \ 

+ 0.1532 + '- 0.0024 , 

170 MeV V 0.1 y' 

which, choosing the central values of the couplings and ai, yields 0.3627, in agreement with 
Tab. 7. 



B Distribution Amplitudes 



In this appendix we collect explicit expressions for some of the twist-3 and 4 DAs that enter 
the LCSRs. These expressions are well-known and have been taken from Ref. [14]. The 
twist-2 DAs have already been discussed in Sec. 4.3. We also motivate and justify the use 
of models for DAs based on a truncated conformal expansion. 

Before defining the DAs, we introduce the hght-like vectors in which they are expressed. 



We denote the meson momentum by (with P 



and the separation between fields 



in a nonlocal operator by (with close to 0) and introduce light-like vectors p and z 
such that 



- P _ 1 H3: 



a;,. 



xP 



{xpy 



x^mr 



/m 



(A.2) 
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m 






m 






yBq-^p 


0.1092 


0.2139 


^0 


0.2036 


0.0990 




0.0867 


0.1552 




0.0467 


0.1743 




0.1034 


0.1641 


rpOq^p 
^3 


0.0303 


0.1455 


yBs^K' 


0.1275 


0.2289 


aBs^K* 

^0 


0.2469 


0.1532 




0.1022 


0.1641 




0.0445 


0.1684 


rpBs^K* 


0.1215 


0.1737 


rpBs^K* 


0.0211 


0.1339 


yBq^K* 

^Bq^K 

B —*K* 


0.1415 
0.1034 
0.1301 


0.2234 
0.1545 
0.1665 


ABq^K 

^0 


0.2071 
0.0614 
0.0436 


0.1269 
0.1658 
0.1386 


yBq^UJ 


0.1048 


0.1884 


^0 


0.1967 


0.0838 




0.0834 


0.1357 




0.0443 


0.1536 


rpBq — >UJ 


0.0984 


0.1440 


rpBq — *ij 


0.0288 


0.1264 


yBs^(t> 


0.1594 


0.2748 


^0 


0.2647 


0.2098 




0.1226 


0.1884 


^2 


0.0558 


0.1784 


rpBs^cp 


0.1469 


0.2019 


rpBs—*(p 

-'a 


0.0316 


0.1433 



Table A: Contributions in and to the form factors at — 0. The numbers correspond 
to the central values of parameter set 2, i.e. = 4.8 GeV. 72(0) follows from Ti(0) = 72(0). 



m 




m 




yBs^K' 
J^Bs^R* 

rpBs-^K* 


-0.0057 -0.0396 
-0.0057 -0.0276 
-0.0056 -0.0297 


aBs-^K* 
^0 

j^Bs^K* 

rpBs^K* 

^3 


-0.0398 0.0024 
0.0079 -0.0394 
0.0104 -0.0293 


yBq^K" 
^Bq^K* 

rpBq^K* 


0.0060 0.0403 
0.0059 0.0281 
0.0059 0.0303 


ABq^K' 
^0 

rpBq^K* 

^3 


0.0403 -0.0001 
-0.0080 0.0395 
-0.0103 0.0299 



Table B: Contributions of Oi to the form factors at = 0. Parameters like in Tab. A. 
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The meson polarization vector e^^'* is decomposed into projections onto the two hght-hke 
vectors and the orthogonal plane as 

We also need the projector onto the directions orthogonal to p and z: 

giv ^Qfiu-— {PtiZu + PvZf,). (A.4) 
pz 

The dual gluon field strength tensor is defined as G^j, = \^f_iupcrG^" . We use the standard 
Bjorken-Drcll convention [42] for the metric tensor and the Dirac matrices; in particular 
75 — ^7^7 ^7^7^, and the Levi-Civita tensor e^j^Ao- is defined as the totally antisymmetric 
tensor with 60123 = 1- This convention differs in sign by the one of Itzykson/Zuber [43] 
used in some programming packages, e.g. FeynCalc. Wc use a sign-convention for the 
strong coupling g where the covariant derivative is defined as — d^ — igA^j^ and hence the 
Feynman-rule for qqg vertices is +ig'Jn. 

Let us also clarify the treatment of SU(3)-breaking effects in DAs. SU(3) breaking occurs 
in three different ways: 

• the contribution of odd Gegenbauer- moments ai^a^... to the DAs of the K*; 

• a difference in the values of the couplings , the even Gegenbauer-moments 03 7^ af^* 
and 3-particle matrix elements; 



the modification of relations between DAs by terms in m^j ± m 



q2- 



We will take into account the first effect wherever it occurs, except for terms in 0(5^), the 
reason being that the structure of 5^ terms is very involved and there are yet unknown 
contributions in rriyaf' induced by 3-particle twist-4 DAs. The second effect is taken into 
account for the decay constants and parametrized by the dependence of the form factors 
on the parameters A, as discussed in Sec. 4.3; wc do not include SU(3) breaking for the 
3-particle matrix elements as information on these effects is virtually nonexistant. The third 
effect is taken into account at 0{Sag,Sas), i.e. for the chiral-even DAs It does not 

occur at 0{5'^) and the corresponding terms are unknown at 0(5^). 

The 2-particle DAs have been defined in Eqs. (15) to (18). Up to twist-4 and 0{5'^), 
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there are seven chiral-odd 3-particle DAs which can be defined as [14] 
{0\q2(z)ac,pgGi,^(vz)qi(-z)\V(p,X)) = 



+ fv^llPae^llgpu - Ppe^±l9L - Pae±lgpi,+ Ppe^±lga^,]Ti^^ 

, rT 2 I (A) ± (A) ± (A) ± , (A) ±irp(4)/ x 

+ fv^v\Pi^(^la9f3u -Pi^e\^9au-Puela9p^^ + Puel^9a^]Ti {v,Pz) 

, fv^V! (A) (A) (A) , (A) 1^(4)/ X 

+ \PaPfie]^pZ„ - ppp,j,e\^z^ - PocPveYpZ^, + ppp^e\^z^,\T^ [v,pz) 

. fv''^V[ (A) (A) (A) , (A) 1^(4)/ V /. r\ 

H \PaPfiel^Zf3 - ppp,j,e\lza - PaPuey zp + ppp^ey Za\Tl [v, pz) , (A.5) 

pz 

{0\q2(z)gGi,^(vz)qi(-z)\V(p,X)) = ifvmyie'ilp^ - e']^lpf,]S(v,pz), 

{0\q2{z)igGf,^{vz)^5qi{-z)\V{p,X)) = ifv'mv[(^±jj>u - e^^lpi^]S{v,pz). (A.6) 

Of these seven amphtudes, T is of twist-3 and the other six of twist-4; higher twist terms 
are suppressed. In the above equations, we use 

T{v,pz) = J Pae-*P^("^-"i+^'^3)r(a), (A.7) 

etc., and a is the set of three momentum fractions a — {ai, a^, a^}- The integration measure 
is defined as i i i 

J Va = J dai J da2 J das 5 ^1 — ^ . (A. 8) 

As for chiral-even DAs, to order 0{5'^) only the twist-3 DAs contribute, which we define as 

{0\q2{z)gGi,^,{vz)'ya'y5qi{-z)\V{p,X)) = fvmvPa\pue''ll - Pf,e';^l]A{v,pz) + 0(my), (A.9) 

{0\q2{z)9G,,{vz)t^^q,{-z)\V{p)) = fymvPalpJH - P^e^llMv,pz) + 0{ml).{A.W) 

At first glance, the sheer number of different DAs, 2 of twist-2, 7 of twist-3 and 9 of twist- 
4, seems to preclude any predictivity of the LCSRs. Appearances are deceiving, though: not 
all these DAs are independent of each other, and one can disentangle their mutual interde- 
pendencies using the QCD equations of motion, which results in integral relations between 
different DAs, e.g. the chiral-odd DAs (f)\\, g^l'^\ gz etc. We shall see examples of such rela- 
tions below. The other important organising principle for DAs is conformal expansion, i.e. a 
partial wave expansion of DAs in terms of contributions of increasing conformal spin. Con- 
formal expansion relies on the fact that massless QCD displays conformal symmetry^"^ which 
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See Ref. [15] for a review on the use of conformal symmetry in QCD. 



37 



allows one to organise the DAs in terms of irreducible representations of the correspond- 
ing symmetry group SL(2,]R). The coefficients of these different partial waves renormalize 
multiplicatively to LO in QCD, but mix at NLO, the reason being that the symmetry is 
anomalous. 

As mentioned above, the plethora of vector meson DAs is not mutually independent, but 
related by the QCD equations of motion. These relations are discussed at length in Ref. [14], 
whose formulas we adapt to the present case. The chiral-even twist-3 DAs are of order 5, so 
we keep the full dependence on terms induced by 0||, but use conformal expansion for the 
admixture of 3-particle DAs: 



(1 - S-,) # = 



u 



u I dv — ^ 



+ u dv 



3_ 
16 



91 



+ 10C3( l-T^c^a"" 



1 
4 



dv 
V 



+ ^u;^^ {5{2u-ir-l}, 

+ 5C3{3(2«- 1)2-1} 



(A.ll) 



+ ^ Ca {^ou^ - (3 - 30(2h - 1)^ + 35{2u - 1)^) , 



(A.12) 



with ^(«) = 20||(77,) + d+{2u - l)0i(M) + d-(l>'Au), 6± = if^/fv) K, ± mgj/my. The 
dimensionless coupling (3 is defined by the (local) matrix element 



{^\(l2gG^^jal5qi\V{P,\)) = fvmvCs 



+ ^fvrrivQ 



= (A) 



(A.13) 



where ^4 is a matrix-element of twist-4. ' ' are matrix elements of quark-quark-gluon 
operators involving derivatives and defined in the second reference of [14]. 

The chiral-odd twist-3 DAs, on the other hand, are 0{5'^), so we model them in conformal 
expansion truncated after the first non-leading order: 



(A.14) 



h\l\u) = 3{2u - if + -a^ {2u - if {5{2u - if - 3) 

+ ^Csc^a'^lS - 30(2ii - If + 35(2ii - 1)^). 
io 



(A.15) 



As mentioned above, we drop contributions in the odd Gegenbauer-moment a^*", as not all 
m^a^ terms are known. 
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IGeV 
2.2 GeV 


0.032 ±0.010 -2.1 ±1.0 3.8 ±1.8 
0.018 ±0.006 -1.7 ±0.9 3.6 ±1.7 



Table C: 3-particle parameters of chiral-even distribution amplitudes. 



/i 


T 

^3 


AT 


/-T 






IGcV 


7.0 ± 7.0 


0.10 ±0.05 


-0.10 ±0.05 


-0.15 ±0.15 





2.2 GeV 


7.2 ± 7.2 


0.06 ± 0.03 


-0.06 ±0.03 


-0.07 ±0.07 






Table D: 3-particle parameters of chiral-odd distribution amplitudes. Terms in are treated 
as described in Sec. 4.3. 



As for the 3-particle twist-3 DAs, we have, quoting from Ref. [14]: 

V(a) = 540C3a;^(Q;i — a-2)a.\a20^i 

1 



l + ^3^_(7a3-3) 



T(a) = 540C3C^3'(tti — Ci2)Ci\Ci2(y\ 



(A.16) 
(A.17) 
(A.18) 



These expressions are valid to NLO in the conformal expansion. 

The chiral-even 2-particle DAs of twist-4, 5(3 and Ay in Eq. (15), are 0(5^), so we drop 
them. For the chiral-odd twist-4 DAs h-^ and Aj^ we use NLO conformal expansion (with 



0): 



hM ^X^\la^-l- 10(CJ ± a)\cT{^u -V)A-la^-^l (a^J \cT{2u - 1) 



(A.19) 



18 



^ + y a - f cr) + (1 4 + ^ ) cTi^u - 1) 



11 2 2 ^ 



126 



70 
11 



(3)N 



X {uu{2 ± 13mm) ± 2^2(10 - 15'u ± 6^^) Inw ± 2#(10 - 15'u ± 6m^) Inw) . (A.20) 

The formulas for chiral-odd 3-particle DAs of twist-4 are rather lengthy and we refrain from 
reproducing them here. They can be found in the second reference of [14]. 

The numerical values of 3-particle matrix-elements are given in Tabs. C and D, for the 
scales 1 GeV and y/rt^^-y^ — 2.2 GeV. The corresponding one-loop anomalous dimensions 
are also given in [14]. The numerical values for the decays constants fy^ are collected in 
Tab. 3. 
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C 3-Particle Contributions to the LCSRs 

In this paper we include contributions of 3-particle DAs to the correlation function (12) at 
tree level. This appendix contains explicit formulas for these contributions. 

The 3-particle DAs of twist-3 have been defined in App. B; the definitions for twist-4 
DAs can be found in Ref. [14]. Their contributions to the correlation functions are most 
easily calculated in the external field method proposed in Ref. [44] . The light-cone b quark 
propagator in an external field reads, in the Fock-Schwinger gauge Xij,A^{x) — 0: 

{0\Tb{x)b{0)\0)A = tSi'\x) + tSi'\x, 0), (A.21) 

with 



Si'\x,0) ^-J ^e-''^-^ j\v^{vS2{k,m)a,,G^%vx)+va,,G^%vx)S2{k,m)^ 



(A.22) 

where Sn{k, m) = {Ij^ + m)/A"'{k) with A(A;) = 1/ (/c^ — m'^)}^ This expression is equivalent 
to Eq. (2.25) in Ref. [4]. The decomposition (29) selects the chiral odd DAs (A. 5), (A. 6) and 
the chiral even DAs (A. 9), (A. 10). Terms in e*^xp/px are treated by partial integration; we 
have checked that all boundary terms vanish. Upon partial integration, we hence have 



px 



I e'^<'-'^f{a,,as)S2{k,m)^ f e'^<''-'^ f{a,, as)Us3(k,m)elkp - el^pA(k)^] 

J a J a 

with l = q + {ai + va3)p, f{x, y) = daf{a, y) and = da^ dai. 

The contribution of 3-particle DAs to the correlation function (12) then reads: 

'-fvmv j^dv j DaA{lf{pq){V{a)+A{a))2vix[TfTf-f^] + 0{ml) 

+ '-fWyQ^v j L'aA(0'5(a)(^;tr[r(g^ + m)^*^^75]+^;tr[r^*j^(^ + m)75]) 
- j^dv j DaA{lfS{a){vix[T{i + m)firi5]+vii[Tfi){-i + m)-i^]) 
^ j^dv j DaA{lfviY[V{4^m)fin^]T^''\o^uOi^) 
+ j^dv j DaA{l)Hv{pq)ir[V{-i + m)fin^\Ti^\oli.Oi^) 
dv j DaA{lfvii[V{i + m)fi>^^]Ti^\a^,a:,) 



^^Note that S^'^^{x,Q) ^ 5(^^(0, —x) as the Fock-Schwinger gauge breaks translational invariance. 
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- dvj DaA{l)Hvipq)tT[rf^{-,i + m)j,]Tt\di,as) 
+ dv J Da {16v{p{e*q)} + 4v{e*pq}-4vm{e*p}){pq)A{lfTl'^\di,a3) 

+ ((v - v){qe*p} - Av{e*(pq)} + v{e*pq} + m(v + 2v){e*p}) A{lfT['^\a^, a^) 
+ dv J Da {-16v{p{e*q)} - 4v{qe*p} - 4:vm{e*p}) ipq)A{lfT^'^\di,a3) 

+ ^((t; + v){qe*p} - Av{e*{pq)} + v{e*pq} + m{v + 2v){e*p}) A{lY T2'^\ai, as) 
In the above formula, we use {abc} — tr[r^^5/75] and {a{bc)} — b ■ c tr[r5i75]. 
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